Experimental investigation on machining characteristics of titanium processed using electrolyte sonicated µ-ECDM system

Micromachining of difficult-to-machine materials is of prime focus nowadays. One such material is Titanium, which has numerous applications in aerospace, chemical, and biomedical industries. The micromachining of Titanium has become the need of the day because of its exhilarating properties. This investigation employs a tailor-made electrolyte sonicated micro-electrochemical discharge machining (ES-µ-ECDM) system to generate microholes in a commercially pure titanium plate with a thickness of 1000 µm. The machining chamber is the ultrasonication unit (36 kHz) with process parameters voltage (V), concentration (wt%), and duty factor (DF) chosen at three levels. The FCC-RSM-based DOE is selected for experimentation to study the machining characteristics like material removal rate, overcut, and circularity. Through holes were achieved at parameters of 80 V, 25 wt%, and 60% DF and 80 V, 30 wt%, and 50% DF. The incorporation of ultrasonication into the system enhanced electrolyte replenishment and evacuation of the debris at the machining vicinity. The assistance technique improved the gas film stabilization around the tool enabling uniform machining. The multi-response optimization is performed using the MOJAYA algorithm to obtain Pareto optimal solutions, and the MADM (R-method) is employed to obtain the optimal parameter. The optimal parameter was found to be 69 V, 30 wt%, and 50% DF, at which the machined microhole was found to have a circularity of 0.9615 with minimal surface defects.

. Conventional µ-ECDM process mechanism for conductive workpiece (a) schematic (b) magnified view of the machining zone (c) formation of the gas bubbles (d) formation of gas film thickness (e) occurance of spark. www.nature.com/scientificreports/ Applying the scale effect, Mouliprashanth and Hariharan 36 carried out electrochemical micromachining. The input parameters voltage, duty ratio, and tool feed rate were varied independently in these studies. Using variations in the pace of material removal, the degree of circularity, and the amount of overcut, they investigated the scale impact at several scales of severity. Similarity theory allowed them to see the scale impact by comparing the granularity of micro and macro systems. In a quantitative analysis of the scale effects, the degree of similarity accuracy was shown to be higher. Anasane and Bhattacharya 37 machined titanium using the ECM process with seven different electrolytes to investigate the anodic dissolution. Among the seven electrolytes chosen five have shown better anodic dissolution. Out of the five three were aqueous electrolytes of the acid, base, and salt whereas the remaining two were non-aqueous base solutions. It was found that the non-aqueous base solution (ethylene glycol with sodium bromide) produced better results and can be considered a cost-effective and ecofriendly electrolyte for machining titanium using ECMM.
Mouliprashanth and Hariharan's 38 research helps to determine and evaluate the geometrical and performance characteristics of ECMM on SS304 alloy under the influence of three distinct kinds of electrolytes: non-passivating, passivating, and composite (CPE). Circularity, conicity MRR, and overcut are geometrical and performance characteristics. The results showed composite electrolyte results were better. The Operandi of optimization with the technique with TOPSIS approach for multiresponse characteristics based on MCDM for ECMM process was also performed to obtain optimal process parameters.
Azhar et al. 39 performed EDM micromachining on three different materials namely aluminum, copper, and stainless steel. They proposed a model (CANFIS) to predict multiple µEDM performances. This model took properties of the material (thermal conductivity, melting point, and electrical resistivity) into consideration to predict MRR, total discharge pulse, overcut, and taperness with more than eighty percent accuracy. Gong et al. 40 performed EDM milling of sintered NdFeB and studied its effect on surface morphology, cross-section, MRR, surface roughness, and electrode wear rate. Pulse current, pulse interval, and pulse duration were taken as process parameters. The EDM milling was compared with EDM and it was observed that EDM milling had better performance on machined surfaces. EDM milling had a large crater radial size, lower surface roughness, thinner recast layer, and less heat-affected zone. The impact of boron carbide B4C powder blended with dielectric fluid on MRR, EWR, and Ra was experimentally investigated by Kolli and Kumar 41 . The addition of B4C powder to the dielectric fluid resulted in a significant increase in the MRR, but a reduced surface finish. Additionally, at high B4C powder concentrations, adhering phenomena were noticed on the machined surface.
Most of the research using the µ-ECDM process focuses on non-conducting materials, and only a few pieces of literature are available for micromachining conductive materials 11 . The following section focuses on the works of literature available on the micromachining of different materials using the µ-ECDM process, particularly conducting materials.
The machining of non-conducting materials has been explored to a vast extent; hence, only a few pieces of literature are discussed in this section, followed by literature corresponding to conducting materials. Torabi and Razfar 42 investigated to examine the ECDM process's potential for the generation of microchannels on polydimethylsiloxane (PDMS). The ECDM technique was capable of fabricating channels on PDMS with greater depth but poor surface quality. Nevertheless, with a reduction in tool diameter, there was a reduction in breadth and surface roughness. Furthermore, grooved tools have produced better outcomes than basic tools. Rathore and Dvivedi 30 used an ECDM process with a sonicated tool on borosilicate glass to further examine the processing capacity. With the tool's sonication, the machining quality was increased, and the thickness of the gas film was better managed. While the MRR rose by 11.13% and penetration depth by 27.17% compared to the standard ECDM process, the reduction in hole overcut was 23.10%.
Oza et al. 43 performed machining on quartz using traveling wire electrochemical discharge machining using zinc coated brass wire. Electrolyte concentration, wire speed, and applied voltage were taken as input parameters. Taguchi's based L9 orthogonal array was used to design the experiments Signal-to-noise ratio and ANOVA were used for the optimization of MRR and kerf width. The findings showed that coated brass wire enhances surface finish and reduces wire breakage. Nema et al. 44 performed machining on quartz material using the ECDM process. A tungsten tool was used to generate micro holes. Feed rate, electrolyte concentration, and voltage were considered as control parameters while objective parameters were maximizing material removal rate and minimizing tool wear. A heat transfer search algorithm was used for optimizing the process parameters to obtain a set of Pareto optimal solutions. Kumar et al. 45 investigated wire ECDM process parameters for machining of quartz material using the Taguchi-based optimization method. According to the findings of the experiments, the zinc-coated brass wire with a 150 µm diameter produced the highest material removal rate and produced the best surface quality. Machining quartz beneath a zinc-coated brass wire may improve surface quality features and the rate at which material is removed. Oza et al. 46 worked on magnetohydrodynamics to improve electrolyte flushing. Microslits were fabricated using magnetic assisted traveling wire ECDM. The magnets generated Lorentz force which applied torque on electrolyte and resulted in a stirring effect. This increased MRR and cutting length by 85.28% and 48.86%, respectively, and reduced the surface roughness by 30.39%.
Huang et al. 47 used the ECDM approach to evaluate the effects of rotating tools (high-speed) on tool electrodes and stainless steel workpieces. Burr-free holes were machined when the TWR (axial) decreased as the tool rotation speed increased. Zhang et al. 48 carried out a comparative study by machining nickel-based superalloy using EDM and Electrochemical Discharge Drilling (ECDD). They investigated the machined holes' surface integrity (surface topography, surface profile, recast layer, surface roughness, etc.). The authors concluded that the surface integrity improved with ECDD with no recast layer and melted debris. Kumar and Das 11 compared the machining of titanium alloy using the EDM and ECSM processes. The MRR and surface integrity of the machined holes were evaluated and found that the MRR and the surface topography were better for ECSM than the EDM process. Also, a thicker recast layer and HAZ were observed using EDM. The crater marks, pit holes, www.nature.com/scientificreports/ and globules were observed to be minimum in the case of ECSM as compared to EDM. Singh et al. 26 employed the µ-ECDM method with a rotating tool mode, metal matrix composites were machined. It was discovered that the L/D ratio was 23.07 times greater than the µ-ECDM without tool rotation. The study found that 900 rpm was the optimum tool rotation for maintaining a consistent gas film thickness. The above literature on ECM and EDM shows the possibility of machining conducting materials with their assistance techniques to improve the machining characteristics. Though µ-ECM and µ-EDM have proven their existence in machining conducting materials the process possesses some merits and demerits which are discussed earlier in this section. The demerits can be addressed by processing the materials using µ-ECDM 12 . The review of the literature using the µ-ECDM process suggests that very few studies have been explored in the area of micromachining of conducting materials, and most of the studies are carried out to understand the surface integrity of the machined surfaces, and the literature lacks an understanding of the machining characteristics such as MRR, OC, circularity, etc. Though the µ-ECDM process has the potential to improve the machining characteristics, it lacks machining beyond 300 µm because of inadequate electrolyte supply at the machining zone. Researchers have suggested some assistance techniques to overcome this complexity. One such assistance is electrolyte ultrasonication, which can replenish electrolyte in the machining zone and enhances the brighter prospects of the µ-ECDM process for machining exotic materials. The main objective of the proposed study is the micromachining of a commercially pure titanium plate of 1000 µm thickness using a tailor-made Electrolyte sonicated µ-ECDM (ES-µ-ECDM) system and optimizing the process parameters for maximizing MRR and circularity and minimizing OC using the MOJAYA algorithm coupled with MADM (R-method).

ES-µ-ECDM prototype
The in-house developed ES-µ-ECDM system is used to carry out the experiments. The 3D schematic view and photographic view of the ES-µ-ECDM system are shown in Fig. 2 and Fig. 3. The system consists of a 36 kHz frequency ultra-sonication chamber used as a machining chamber with fixtures assembled to hold the workpiece. A regulated DC power supply (0-120 V and 0-5 A) is employed in the system with a pulse generation circuit assembled along with it. The pulse generation circuit helps in producing pulsating DC of the required frequency and duty factor (0-100%). A digital storage oscilloscope is connected to the pulse generation circuit to understand whether the required frequency and duty factor are attained or not. The developed system also consists of a three-axis motion machine driven by stepper motors connected to stepper motor drivers, controlled by a universal g-code sender (UGS) coded using Arduino UNO. The SMPS power sources provide 12 V constant DC to stepper motors, Arduino UNO, and exhaust fan. A worm wheel and worm arrangement are also aligned to the axis of the stepper motor to reduce the feed rate.

Experimental procedure
A detailed experimental study is conducted to understand the effects of the process parameters on the machining characteristics. In-depth pilot experiments are performed to finalize machining process conditions to accomplish the present task. A commercially pure titanium workpiece of 1000 µm thickness is machined with a titanium tool of 800 µm diameter in NaOH solution. The NaOH electrolyte has a conductivity of 414 mS/cm at a concentration of 20 wt% which reduces to 292 mS/cm with the increase in concentration to 30 wt% at 25 °C. Though at room temperature the conductivity is less during the ECDM process the conductivity increases with rising in temperature as it is known that the kinetic energy of the particles increases 49 . The tool and workpiece are cleaned using acetone to remove the impurities and dirt on the surfaces. The tool is connected to the negative Pilot experiments. The experiment was conducted initially with 20 wt% concentration and 50% DF, and the voltage was kept increasing till the spark occurred, and it was observed that the spark took place at a voltage of 35-37 V. The voltage was set to 40 V, and machining was conducted. It was observed that there was no material removal, and the tool bend when it came in contact with the workpiece showed only a heat spot on the surface of the workpiece. The voltage was then increased until an appreciable material removal was achieved, which was found to occur at 60 V. Keeping the concentration and voltage constant at 20 wt% and 60 V and reducing the DF to 40% DF led to material removal occurring along with tool bending, which is not desirable as it affects the circularity and taperness of the machined microhole. As a result, the lower level of the process parameters is fixed to 60 V, 20 wt%, and 50% DF. Then the concentration was increased to 30 wt%, and the DF was kept constant at 50% DF. The voltage was also increased to 70 V and 80 V, which enabled higher material removal and higher overcut. The increase in voltage above 80 V resulted in lower MRR because of the melting of the tool (axial wear). The higher voltage level was fixed to 80 V; now, the concentration was increased to 35 wt%, which resulted in higher electrolyte fumes with higher TWR, which is not desirable. The higher level of concentration and voltage was set to 30 wt% with 80 V. Now, keeping concentration and voltage constant at 30 wt% and 80 V, the DF was increased to 60% DF and then 70% DF. At 70% DF, there was enough material removal with acceptable TWR, and above 70% DF, higher axial tool wear reduced the MRR. Hence the higher level of the factor was set at 80 V, 30 wt%, and 70% DF.
From the pilot experiment and the extensive literature review, the variable and constant process parameters for a detailed experimental study were finalized and tabulated in Table 1.
A face-centered central composite RSM (FCC-RSM) design of experiments (DOE) was chosen to carry out the experiments. A three-level three-factor experimental setting has 3 3 = 27 experiments. The full factorial DOE of FCC-RSM enables the experimental layout with 20 experimental runs (14 axial points and six central points). The machining characteristics are MRR, overcut (OC), and circularity. These are used in the performance evaluation of the process. A multi-objective metaheuristic optimization was carried out using the MOJAYA algorithm paired with the MADM (R-method) to acquire the best process parameters. This optimization aimed to maximize MRR and circularity while reducing OC. The MRR was estimated using a weighing balance (Make: Mettler Toledo) by measuring the weights before and after the machining using Eq. (1). The OC and circularity were estimated using an optical microscope with image processing software (ZEISS Stemi 2000-C). The circularity at the entrance of the machined holes is measured as the ratio of minimum to maximum Feret diameter. Figure 4 illustrates the methodology involved in calculating the circularity at the entrance. The equation Eq. (2) and Eq. (3) is used to where D min and D max are the minimum and maximum Feret diameters of machined holes at the entrance.

Results and discussion
The machining characteristics corresponding to the FCC-RSM layout are tabulated in Table 2. The machining was conducted two times, and the average values were considered.
Material removal mechanism in ES-µ-ECDM system. The detailed material removal mechanism involved in the ES-µ-ECDM system is shown in Fig. 5. The ultrasonication unit is used as an electrolyte chamber where the base plate of the ultrasonicator is attached to a transducer. The transducer provides very subtle vibrations when connected to the power supply. The motion of the base plate has two strokes. The downward stroke (rarefaction) of the base plate creates lowpressure zones inside the electrolyte fluid. These low-pressure zones initiate the boiling of electrolyte resulting in micro-size gas bubbles (voids) because of the cavitation phenomenon as shown in Fig. 5f 50 . The voids created due to the downward stroke of the base plate are of microscopic scale and cannot be seen during operation. During the upward stroke (compression), the pressure in those regions increases and reduces the size of the gas bubbles. After passing through several cycles of compression and rarefaction the bubbles tend to grow in size reaching  www.nature.com/scientificreports/ an unstable state as shown in Fig. 5e. At this state on the next compression stroke, the bubble implodes at a high localized jet velocity 51 . This phenomenon is continued in the electrolyte during the whole process but its effect is visualized during pulse off time of the machining. The discharge action takes place during the pulse on time and soon after discharge there becomes a zone of vacuum with very low pressure. During this time, the constant implosion of the voids near the machining zone helps in flushing the electrolyte with high velocity because of the higher pressure difference created. This enables the removal of debris from the machining vicinity and also produces fresh electrolyte for the process to continue. Figure 5a-d explains the procedure of generation of the spark which is quite similar to the conventional µ-ECDM process. In the conventional µ-ECDM process after the hydrodynamic region is reached due to the lack of the electrolyte only side sparking occurs which increases the entrance hole overcut but does not aid in the depth of penetration and also protrusions are formed at the  Influence of voltage on machining characteristics. The ES-µ-ECDM system is used for the sonication of the electrolyte, and voltage plays a prominent role in the micromachining of commercially pure titanium workpieces. Stabilizing the spark and spark intensity depends on the gas film's thickness and voltage, aiding assistances like sonication 52 . The effect of voltage on the machining characteristics is shown in Fig. 6. It can be observed that increasing the voltage from 60 to 80 V increased all the machining characteristics. Ultrasonication provided to the electrolyte helped reduce the gas bubble's departure radius, resulting in gas bubbles with a thin and more stable gas layer. The latter leads to the channelized spark under the tooltip and thus, increasing the depth of penetration and generating through-holes at higher levels of voltage. The rise in voltage increased the spark intensity at the machining zone enabling a higher MRR. As a result of the increase in voltage from 60 to 80 V, an average percentage increase in MRR of 89.67%. The increase in spark intensity also increases the OC and circularity associated with the machined hole. The percentage increase in OC and circularity with an increase in voltage from 60 to 80 V, keeping concentration and duty factor constant at 20 wt% and 50% DF, is 73.91% and 0.44%, respectively. The average percentage increase in OC and circularity with a surge in voltage is 64.09% and 3.71%, respectively. The SEM images in Fig. 7 show the rise in OC and circularity. As the spark intensity increases, there is an evident increase in OC, but the increase in hole circularity might be due to proper melting and vaporization of the workpiece at higher voltages. The rise in OC also states the evidence for a rise in MRR. It is also observed that with an increase in voltage the HAZ increases but the formation of the passivation layer was seen to be less. This might be because higher spark intensity due to a surge in voltage increased the HAZ and does not allow the formation of a passivation layer. As the spark hindered the formation of the passivation layer the depth of penetration increased and through-holes are achieved at a higher voltage of 80 V.
Influence of electrolyte concentration on machining characteristics. The sonication of electrolytes helps flush debris out and circulate fresh electrolytes at the machining zone. The increase in electrolyte concentration increased the number of conducting ions, thus, increasing the spark intensity. The aiding of sonication to the system increased penetration depth, which eventually resulted in through-hole formation. As shown in Fig. 8, the plots elucidate the effects of electrolyte concentration on the machining characteristics. The increase in concentration from 20 to 30 wt% enhanced all machining attributes. The surge in MRR, OC, and circularity can be attributed to the increase in anodic dissolution and intensity of the spark because of an increase www.nature.com/scientificreports/ in the number of conducting ions in the machining zones vicinity. The rise in OC and circularity can be justified by the SEM images shown in Fig. 9. It can be observed from Fig. 9a and b, that with a surge in concentration from 20 to 30 wt%, keeping voltage and duty factor constant at 70 V and 50% DF, the percent increase in OC and circularity is found to be 3.13%, and 3.71%, respectively. The average percentage surge in MRR, OC, and circularity with increased concentration is 42.45%, 12.52%, and 3.73%, respectively.
Influence of duty factor on machining characteristics. The duty factor is one of the critical parameters in machining, as proper tuning helps produce better-machined surfaces. Continuous DC lacks the ability  www.nature.com/scientificreports/ of electrolyte replenishment in the machining zone, enabling increased thermal damage and heat-affected zones on the workpiece 53 . At low machining voltages, pulsing DC is responsible for making smooth surfaces. The plots in Fig. 10 indicate the influence of the duty factor on the machining attributes. The higher duty factor boosted MRR and OC while decreasing circularity. The rise in MRR and OC might be attributed to the existence of a discharge for a longer time, resulting in greater material removal at the surface leading to higher OC. The reduction in circularity of the machined microholes might be attributed to a longer pulse on time, which does not provide enough time for the molten material to cool, causing it to adhere to the machined hole's surfaces and resulting in irregular micro holes. The average percentage surge in MRR and OC due to a rise in duty factor is 2.10% and 4.27%, respectively. The average percent reduction in circularity is found to be 3.59%. This percentage increase in  www.nature.com/scientificreports/ OC and reduction in circularity can be justified by the SEM images shown in Fig. 11, and it can be observed that with an surge in duty factor from 50% DF to 70% DF keeping voltage and concentration constant at 70 V and 25 wt%, the percentage rise in OC is 13.67%, and the percentage reduction in circularity is found to be 4.14%. The increase in duty factor also increased the HAZ as the duration of the discharge is more. From the above discussion, it can be established that voltage is the major factor, followed by concentration for the machining characteristics like MRR and OC. Nevertheless, the circularity of the machined microholes is influenced by all the parameters with an edge-to-duty factor. The latter influences the circularity and is more effective than voltage and concentration. The SEM image shown in Fig. 12 elucidate that at 80 V, 30 wt%, and 50% DF through-hole was obtained, but at 80 V, 30 wt%, and 70% DF, a through-hole was not obtained, which might be due to a higher tool wear rate (axial wear). It is also found that the circularity of the machined hole decreased with an increase in duty factor. The circularity of the hole machined at 50% DF, keeping voltage and concentration at 80 V and 30 wt%, is higher than at 70% DF. The ultrasonication of the electrolyte helped in the proper replenishment of electrolyte from the machining zone and the successful evacuation of the debris formed. The channelization of the spark enabled uniform machining and the SEM micrographs shown in this article clearly show that no such protrusions are formed.
Mathematical modeling and optimization. For the regression study, mathematical modeling was done using Design Experts-13 statistical software. The RSM modeling found quadratic relationships between the process parameters and the machining characteristics. Equations (4)-(6) are the respective regression equations for MRR, OC, and Circularity. The ANOVA results of the machining characteristics are given in Table 3, and it can be observed that voltage and concentration are significant for all the machining characteristics. However, the duty factor is significant for circularity and not MRR and OC. This statement agrees with the statement given at the end of the results and discussion. The ANOVA table (Table 3) suggests model adequacy for all the machining characteristics as the model is significant and lack of fit is not significant, which is desirable. The R 2 and Adj. R 2 value for MRR, OC, and circularity are 0.9815, 0.9648, 0.9790, and 0.9601, 0.9916, and 0.9840, respectively. These values are above 95%, which supports the model adequacy.  www.nature.com/scientificreports/

Regression equation.
where A, B, and C represent voltage, concentration, and duty factor, respectively.

MOJAYA Optimization coupled with R-method (MADM). The regression equation obtained from
the RSM modeling is used as the objective function to perform optimization using the JAYA algorithm. The criterion selected for optimization was maximization of MRR and circularity and minimization of OC. A MAT-LAB code was developed for performing the MOJAYA optimization and validated with Rao et al. 54 . The Pareto optimum front was obtained using the posterior version of the optimization. As illustrated in Fig. 13, the Pareto optimal front for MRR, OC, and circularity objectives. A population size of 100 and an iteration or termination criteria of 100 were the parameter values utilized to generate 70 Pareto optimal solutions. MOJAYA algorithm does not have any algorithm-specific parameters that have to be tuned to reach a global optimum. This advantage of MOJAYA helps the algorithm find a global optimum instead of ending up in a local optimum 55 . The 70 Pareto optimal solutions obtained are all optimal as they are non-dominated solutions. To find an optimal value, different MADM techniques can be used, such as VIKOR 56 , TOPSIS 57 , AHP 58 , GRA 59 , PROMETHEE 60 and R-method 61 . This article uses R-method to find optimal parameter among the generated Pareto optimal solution. The responses are ranked based on the criterion of maximizing and minimizing. If the criteria are maximization, then the responses are ranked in decreasing order, and if the criteria are minimized, the responses are ranked in increasing order. The objective ranking is solely dependent on the decision-maker and which parameter  www.nature.com/scientificreports/ is more useful for a particular application. The circularity is considered a major response, ranked 1, trailed by OC and MRR, ranked 2 and 3, respectively. A MATLAB code is developed for R-method and validated by Rao et al. 61 . Equation (7) is used to estimate the weights of the objectives and responses. Once weights are obtained for both objectives and responses, the composite score is calculated by summating the product of weights of objectives and responses. Finally, the row corresponding to the highest composite score is ranked one, followed by the others in decreasing order. Table 4 shows the Pareto optimal solutions, responses, weight calculations, composite score, and rank. The process parameters for rank "1" are considered to be optimal. It can be observed that the optimal process parameters obtained by the integrated MOJAYA algorithm and R-method are 69.36779 V, 29.9294 wt%, and 50% DF. The predicted values are rounded to the nearest integer, and validation tests are conducted to check the feasibility of the results obtained. The process parameters are rounded to 69 V, 30 wt%, and 50% DF, and experiments were conducted five times. The average values of the machining characteristics are tabulated in Table 5. It can be seen that the predicted and the experimental values are in close agreement with each other, with an error percent not exceeding ± 5. The SEM image obtained for the experiment conducted at the optimal parameter setting is shown in Fig. 14. The SEM micrograph suggests that the machined microholes obtained at optimal process parameters are smooth with high circularity and minimal OC.

Conclusions
In this article, the experimentation was conducted on a commercially pure titanium plate of thickness 1000 µm using a titanium tool of diameter 800 µm for machining microholes using a tailor-made ES-µ-ECDM system characterizing the machining responses like MRR, OC, and circularity. Further, multi-objective optimization is accomplished using the MOJAYA algorithm, and the optimal process parameter value is predicted using R-method (MADM). The following conclusion is drawn from this study: • The ultrasonication of the electrolyte helped stabilize the spark by improving the machining depth with an increase in process parameters. Through holes were obtained at two parameter settings of 80 V, 25 wt%, and 60% DF and 80 V, 30 wt%, and 50% DF. • The ultrasonication of the electrolyte enabled uniform machining by channelizing the spark at the machining zone even at higher depths (> 300 µm) and produced protrusion-free surfaces. • The assistance of ultrasonication helped in proper flushing of the debris from the machining zone and replenishes the electrolyte which established a stable gas film around the tool surface which enhanced the depth of penetration.  www.nature.com/scientificreports/ • The rise in voltage and concentration increased all machining characteristics, but with a rise in duty factor, the MRR and OC increased, and the circularity reduced. • The ANOVA of the machining characteristics suggested that the duty factor helps in improving surface characteristics like circularity. • The 70 Pareto-optimal solutions were acquired by performing the posterior version of metaheuristic optimization. The optimal process parameter was obtained using R-method, where circularity was considered the prime objective, followed by OC and MRR. • The optimal value predicted was 69.36779 V, 29.9294 wt%, 50% DF, which was rounded to 69 V, 30 wt%, and 50% DF to perform the experiment, and the results obtained at these parameter settings are in good agreement with the error percent not exceeding ± 5%.  Figure 14. Optical microscope and SEM micrographs of the machined multiple microholes at the optimal process parameter setting of 69 V, 30 wt%, and 50% DF. www.nature.com/scientificreports/